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Abstract: In this paper, we address the aggregation of dependent stop loss reinsurance risks where the de¬ 
pendence among the ceding insurer (s) risks is governed by the Sarmanov distribution and each individual risk 
belongs to the class of Erlang mixtures. We investigate the effects of the ceding insurer(s) risk dependencies on 
the reinsurer risk profile by deriving a closed formula for the distribution function of the aggregated stop loss 
reinsurance risk. Furthermore, diversification effects from aggregating reinsurance risks are examined by deriv¬ 
ing a closed expression for the risk capital needed for the whole portfolio of the reinsurer and also the allocated 
risk capital for each business unit under the TVaR capital allocation principle. Moreover, given the risk capital 
that the reinsurer holds, we express the default probability of the reinsurer analytically. In case the reinsurer 
is in default, we determine analytical expressions for the amount of the aggregate reinsured unpaid losses and 
the unpaid losses of each reinsured line of business of the ceding insurer (s). These results are illustrated by 
numerical examples. 

Key words: Risk aggregation; Sarmanov distribution; Mixed Erlang distribution; Capital allocation; Stop loss 
reinsurance; Reinsurance default risk; Default Probability. 

1 Introduction 

Reinsurance companies operate in many regions in the world and insure various insurance business lines. In 
this respect, it is well recognised that the ceding insurer(s) losses are dependent. This risk dependency can be 
seen between individual risks within each insurance portfolio and also across business lines. Furthermore, the 
phenomena of dependence also occurs from global risk factors which generate claims simultaneously to each 
business line, for instance an hurricane damages buildings or cars which affect property lines, at the same time, 
causes people injuries which influence accident lines. In the risk management framework, for instance the Swiss 
Solvency Test (SST), similarly to insurance companies, reinsurance companies are obliged to hold a certain level 
of risk capital in order to be protected from unexpected large losses. The determination of this capital requires 
the aggregation of the losses generated from each reinsurance portfolio whose distribution depends on the loss 
distribution of the ceding insurer (s). Meyers et al. [14] is one of the first contribution which have addressed 
the aggregation of dependent reinsurance risks to evaluate risk capital. In this regard, in order to derive 
explicit formula for the measure of risk capital including Value-at-Risk (VaR), Tail Value-at-Risk (TVaR) for 
the aggregated risk, an important task is the appropriate choice of the marginals and the dependence structure 
between risks. For our framework, mixed Erlang distribution has been chosen as a claim size model for the 
individual risk of the ceding insurer (s). One of the reason of the tractability of this distribution is the fact that 
the convolution of such risks belongs again to class of Erlang mixtures, see [5]. Thus stop loss and excess of 
loss premiums have a closed expression which are very usefull in reinsurance risk modelling, see Lee and Lin 
[nun]. In this contribution, we address the dependence strncture between risks by the Sarmanov distribution. 
The aim of this paper is to analyse the effects of the ceding insurer(s) risk dependencies on the reinsurer risk 
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profile which has only stop loss reinsurance portfolios. Diversification effects from aggregating reinsurance risks 
are examined by deriving a closed expression for the risk capital needed for the whole portfolio and also the 
allocated risk capital for each business unit. The effects of the reinsurer default are also analysed. The paper is 
organised as follows: in Section 2 we describe the background of the Sarmanov distribution as a model for the 
dependence structure between insurance risk and the mixed Erlang distribution with a common scale parameter 
as a claim size model. The risk model of the ceding insurer is explored in Section 3, with numerical examples, by 
deriving the joint tail probability of the aggregated risk of two portfolios. In Section 4, the aggregation of stop 
loss mixed Erlang risks of a reinsurer is addressed by determining a closed form for the distribution function 
(df) of the aggregated risk. Capital allocation and diversification effects are also presented with numerical 
studies. We also analyse the default risk of the reinsurer by deriving an analytical form for the expected unpaid 
losses and the default probability with numerical illustrations. All the proofs are relegated to Section 5. Some 
properties of the mixed Erlang distribution are presented in the Appendix. 

2 Preliminaries 

2.1 Sarmanov distribution 

Due to its flexibility to model the dependence structure between random variables (rv), the Sarmanov distribu¬ 
tion, introduced in Sarmanov have been widely used in many fields. Concerning insurance applications, to 
calulate Bayes premiums in collective risk model Hernandez et al. [8] have addressed the dependence between 
risk profiles using multivariate Sarmanov distribution. Sarabia and Gomez m have used Sarmanov distribution 
to fit multivariate insurance count data with Poisson-Beta marginals. The contributions [271126] have explored 
tractable asymptotic formulas in the context of ruin probabilities where the dependence between insurance risks 
is governed by the Saramanov distribution. Refering to [10] . a random vector (Wi,..., A„) has multivariate 
Sarmanov distribution with joint density given by 

n / n h \ 

+ X! := (a::i,...,x„), (2.1) 

i=l ^ /i=2 l^ji<:/2<...<lfc^n fc=l 

where 4>i are kernel functions, which are assumed to be bounded and non-constant such that 

E{<(.,(W)} = 0, 

n h 

1 -i- 'y ( y ( 4'jki^jk) ^ Oj G K (2.2) 

h=2 lt^ji<j 2 <...<jh^n k=l 

are fulfilled. Some general methods for finding the kernel function (pi was specified by Lee [10] for different 
types of marginals. In particular, it is commonly used to choose (pi{xi) = gi{xi) — E{gi(Ai)} for marginal 
distributions with support in K-|_ (see e.g. US]). The following three cases are the usual specifications of gi{xi): 

(i) gi{xi) = 2Fi[xi) which corresponds to the Farlie-Gumbel-Morgenstern (FGM) distribution, where Fi is the 
survival function of Aj, 

(ii) gi{xi) = xl — E{A*} such that the t-th moment E{A*} of Xi is finite, 

{Hi) gi{xi) = — E where E < oo is the Laplace transform of Xi at t. 

2.2 Mixed Erlang Marginals 

These last decades, mixed Erlang distribution with a common scale parameter is one of the most usefull model 
for insurance losses. In risk theory, using the mixed Erlang distribution as a claim size model, an analytical form 
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for the finite time ruin probability has been derived by Dickson and Willmot [B] and Dickson . Recently, using 
the EM algorithm, mixed Erlang distribution has been fitted to catastrophic loss data in the United States by 
Lee and Lin m and also to censored and truncated data by Verbelen et al. m . Moreover, Lee and Lin |12) . 
Willmot and Woo |25j have developed the multivariate mixed Erlang distribution to overcome some drawbacks 
of the copula approach while Badescu et al. [I] have used multivariate mixed Poisson distribution with mixed 
Erlang claim sizes to model operational risks. Eurthermore, Hashorva and Ratovomirija [7] have addressed risk 
aggregation and capital allocation with mixed Erlang marginals and Sarmanov distribution. In the sequel, we 
denote respectively 


Wk{x,l3) = 


k ^k—1 


P^x' 


Wk{x,P)=J2 

j=k 


{Pxpe 

j! 


k-l 


Wk{.x,p) = Y^ 


{Pxpe 

j! 


-,x > 0, 


(2.3) 


(fc-1)! 

the pdf, the df and the survival function of an Erlang distribution where fc € N* is the shape parameter and 
/3 > 0 is the scale parameter. As its name indicates, the mixed Erlang distribution is elaborated from the Erlang 
distribution, its pdf and df are respectively dehned as 


f{x,P,Q) = ^qkWk{x,P), F{x,P,Q) = '^qkWkix,P), (2.4) 

fe=i fc=i 

where Q = (qi, q 2 , ■ ■ ■) is a vector of non-negative weights satisfying = 1- Hereafter we write X ^ 

ME{P, Q) if X has pdf and df given by (12.41) . One of the main advantages of the mixed Erlang distribution in 
insurance risk modeling is the fact that many useful risk related quantities, such as moments and mean excess 
function, have explicit expressions, see e.g., [n m mui]. Eurthermore, the mixed Erlang distribution is 
a tractable marginal distribution for the Sarmanov distribution. Next we present a result for the correlated 
insurance portfolios. 


3 Ceding insurance risk model 


In this section, we consider two insurance portfolios which both of them consists of k risks and we denote 
Si^k = Sfci ^ 2 ,k = X)i=fc+i the aggregated risk of each portfolio where Xi,i = I,..., 2fc is a positive 

continuous random variable (rv) with finite mean. Hereafter, we assume Xi ~ ME{Qi, Pi),i = ... ,2k and the 

dependence structure between risks within and across the portfolio is governed by the Sarmanov distribution 
with kernel function 

PpXi) = gPxi) -E(5i(X*)), 

which shall be abbreviated as 

(Ai,...,A2fe)-5ME(/3,Q) 

where (3 = {Pi,..., P 2 k), Q = {Qi, ■ ■ ■, Q 2 k)- In the rest of the paper we consider for gi one of the three cases 
described in (i), {ii) and (Hi). 

Eurthermore, we define two vectors of mixing weights Q{Qi) and 4f(Qi) where their components depend on the 
kernel function pi. In particular, the components of Q{Qi) = {di,i,di, 2 , ■ ■ •) are given by: 


• for gi{xi) = 2F{xi), = 5 ^ Z)j=i 

• for gi{xi) = x\. 


's-V 
kJ - h 


s — 1, 2, ..., 


0 

r(s) 

Y^oo ~ r(j+t) 

2^i = l T{j) 


for s ^t, 
for s > t, 
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• for gi{xi) = e with /3 = ;gf^,s = 1 , 2 ,..., 

2^j=i 91JP ^ 

whilst the components of ^(Qi) = (V'i,!,'0i,2, • ■ •) are given by 


'4’i,s — 'y ] Qi,j 


'k-r 


j=i 


, j — 1/ \Z{f32k) 


|3^ 


1 - 


Z{P2k) 


where Z{l32k) = 2 ^ 2 fc for gi{xi) = 2F{xi), Z{P 2 k) = P 2 k for gi{xi) = x\ and Z{l32k) = P 2 k + t for gi{xi) = e 
Moreover, for given mixing weights Vi = {vii,Vi 2 , ■ ■ ■),% = 1,..., n + 1 we define a vector of mixing probability 
n(lfi, • ■ ■, fon+i) as follows 




0 for I = 1 ,..., n. 


for l = n + l,n + 2 ,.... 

We present next the main result of this section. 

Proposition 3.1 If {Xi, . .., X 2 fc) SME{f3, Q) with 132k > I3i,i = 1,..., 2fc — 1, then the joint tail probability 
of Si^k and S 2 ,k is given by (set ■=^{9jmiZ^jm)}) 

2k i 

F {Suk > ui, S2,k > U 2 ) = fiFs,,k(.ui)Fs^^^{u2) + Y^ Y[ljmFs^ ,^{ui)Fg^_^{u2), 


l■=l 3i,32,...,jl 


m—1 


where 


2k 


Cl “ 1 J2 7ii7i2 ^iij2J3 7ii7i2 7i3 + ■•■ + ( 1) Ji 


31 32 


31 32 J3 


2=1 


Cil — ^iij27i2 + '^ji,i2j37^2 7^3 + ■•■ + ( 1) 72 ; ? 

i2 33 


31 32 


^h,h,33lh fo ,3ilhl3i + ■ • ■ + (-l)^^Q!y,,,,2fe h ), 

31 32 ^ ji 33 3i j6C\{ii,j2} 

(,31,32,33 ~ E! E! E! ( rr ji ,32,33 ~ E! rr ji, j2 ,33,34 7j4 + E! E! rr ji, j2 ,33,34, ,33 Iji 7j5 + ■ . ■ + (~1) CXl,...,2k 7i 

fl J2 33 34 34 fs *GC\{jl J'2.j3} 

0l,...,j2fc-l “ E^ ■ ■ ■ ^ ^ ~ a^l,...,2k'yj2k^ 

31 32k-l 

(jl,■■■,32,1c ~ oil,...,2k, 

with C = {1,.. .,2k}, ji e C,j 2 e C'\{ji},j 3 G C\{ji,j 2 }, ■ ■ ■ ,j 2 k G C\{ji,... ,j 2 fe-i}, 

51, k MEinmQi ),..., ^{Qk), Z(p2k)}), 

52, k ME(IimQk+i),..., mQ2k), Zif32k)}), 

51, k MEiurngi), ■■■, mQi),z{fi2k)}), 

52, k ME{umQU ,),..., ^iQ;,),Z{P2k)}), 


and for i = 1,... ,2k 


q: = 


Qi */ ii{ji,j2,...,ji}, 

QiQi) if i & {ji,j2,...,ji}. 


Example 3.2 Assume that the ceding insurer has two portfolios say Portfolio A and Portfolio B. Concerning 
the dependence structure betwen risks, two cases of kernel function are considered 4>i(xi) = 2Fi(xi) — 1 which 
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defines the FGM distribution as explored in and (j)i{xi) = e — E{e introduced by Hashorva and 
Ratovomirija ^ for mixed Erlang marginals. In the rest of the paper we refer to the latter as the Laplace ease. 
Table [3R\ vresents the parameters of each individual risk X^fi = 1,..., 4 and their central moments, whilst Table 
displays the dependence parameters between Xi, X 2 , X^ and X 4 . We note that these dependence parameters 
have been chosen so that ( 1 ^ holds. 



X, 

A 


Mean 

Variance 

Skewness 

Kurtosis 

Portfolio A 

Xi 

0.12 

( 0 . 4 , 0 . 6 ) 

13.33 

127.78 

1.55 

6.50 

X 2 

0.14 

(0.3,0.7) 

12.14 

97.45 

1.49 

4.33 

Portfolio B 

A 3 

0.15 

(0.5,0.5) 

10.00 

77.78 

1.62 

6.80 

X 4 

0.16 

( 0 . 8 , 0 . 2 ) 

7.50 

53.13 

1.88 

8.16 


Table 3.1: Parameters and central moments of Xi, i = 1, 2,3,4. 



Oil,2 

0^1,3 

0:1,4 

02.3 

02,4 

03,4 

Ol,2,3 

Oil,2,A 

Ol,3.4 

02,3,4 

01,2,3,4 

FGM 

0.6 

0.1 

0.1 

0.1 

0.04 

0.5 

0.11 

0.12 

0.10 

0.15 

0.07 

Laplace 

16 

5 

3 

5 

3 

8 

56 

30 

15 

20 

170 


Table 3.2: Dependence parameters of (Xi, X 2 ,X 3 , X 4 ). 


It can be seen from Table UTSl that the interdependence between the two insurance portfolios yields high probability 
for the aggregated risk of each portfolio to exceed simultaneously some threshold. 


Thresholds 

Independence case 

Laplace case 

FGM case 

iui,U 2 ) 

P (5i,2 > Ml, 82,2 > M2) 

IP (5i,2 > Ml, S 2.2 > M2) 

IP(5i,2 > Ml, 52,2 > M2) 

(20,15) 

0.1494 

0.1569 

0.1573 

(25,20) 

0.0697 

0.0751 

0.0795 

(30,25) 

0 . 0.304 

0.0331 

0.0374 

(35,30) 

0.0125 

0.0138 

0.0165 


Table 3.3: Joint tail probability of 51,2 = Xi + X 2 and 82,2 = X 3 + X 4 . 


4 Reinsurance risk model 

In this section, we denote R 2 := Ti,fc+T 2 ,fc the aggregate reinsurance stop loss risk where Ti,fc := —di)_|_ and 

T 2 ,k '■= (52,fc — ^ 2 )+ represent two stop loss reinsurance portfolios with = Si=fe+i 

the ceding insurer aggregated risk and di,i = 1,2 some positive deductible. Additionally, for a given risk 
X ^ ME{fi,Q) with df F and for a deductible d > 0 we denote in the rest of the paper 

00 

Fx{y + d) = ^ Afe(d,/3,Q)lTfe+i(y,/3), 

fc—0 
00 

Ux{c,d,l3) = '^{k + l)Ak{d, fi,Q)Wk+ 2 {c, 13), 

k —0 


Ak{d,fi,Q) 


1 , ^ 

■n'^<l3 + k+lWj+l{d,l3). 

^ 3=0 


with 
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Furthermore, for Xi ~ ME{Qi, P), with di > 0, z = 1, 2 we define 


Fxi+X2{di,d2, s) = EE Afc(di, /3, Qi)Aj{d2, fi, Q2)Wk+j+2is, /3), 

k—0 j—0 
^ oo oo 

C/xi(c,di,d2,/3) = + l)Afc(di,/3,gi)Aj(d2,/3,Q2)fFfc+j+3(c,/3), 

^ fc-0 j^o 

^ oo oo 

Ux2{c,di,d2,P) = -g E E^'^ gi)Aj (^2, /3, g2)hF fc+j_i_3(c, /?), 

fc— 0 j—O 
^ oo oo 

C/X 1 +X 2 (c, di, d2, /?) = ^ X! E(* + -^ + 2)Afc(di, /3, Qi)Aj{d2,/3, Q2)Wk+j+3{c, /3). 

A:—Oi—0 


4.1 Aggregation of reinsurance stop loss risks 

In the following result we show that the df of the aggregated stop loss risk R2 has a closed form which allows 
us to derive analytical formula of its mean excess function. 


Proposition 4.1 If {Xi,..., X2k) ~ SME{f3, Q) with P2k > ldi,i = 1,..., 2fc — 1 and dj > 0, j = 1, 2, then the 
df of the aggregated stop loss risk R2 is given by 


Fr -2 (s) 


Fs,,^,S2,Adi,d2) for s = 0, 

Fsi,k,S2.kidi+s,d2 + s) for s > 0, 


where 


Fsi,k,S2,k{d^^d2) = Cit^Si,fc(di)Fs2 Jd2) 


2 k 

E 

l=i 


E 0 




(4.1) 


Fsi_k,S2.kidi + s,d2 + s) = + s) + Fs^i^idi + s)Fs^,^{d2) + Fs^,^+S2,kidi,d2,s)^ 


2 k I 

+ E E ^oi,02,...,jiY[ljm[Fs,^{di)Fg^^{d2 + s) 

* = 1 tl Ji 


m—1 




with Si,k, S2,k, Si^k, S2,k are defined in Provosition \3 . 1\ 

Remarks f.2 Given the tractable form of the df in mi, many risk related quantities for R2 have an explicit 
form, for instance, for c > 0 the mean excess function of R2 is given by 

E {i ?2 — c|i ?2 > c} = 


Fr2 (c) 


Cl ( {dl)US2,k (c, d2,Z{fl2k)) + -Fsa.fc id2)U (c, di,Z{f]2k)) 


2 k 


+Us^,k+S2.kic,di,d2,Z{f32k))] +J2 E Yi ljm[Fs,^^{di)U §^^^{c,d 2 , Z{l32k)) 

l=i jl,j2,---,jl m=l 


+F§2 ^ (d2)Ug^^^ (c, di,Z(^2k)) + ^ (c, di,d2,Z{P2k)) 


(4.2) 


Example 4.3 In this illustration, we consider the same parameters of each individual risk of the ceding insurer 
portfolios as in Table \3.1\ Furthermore, we assume that the ceding insurer re-insures its two portfolios to a 
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reinsurer with stop loss programs where the deductibles are di = 40 and c ?2 = 30 for Portfolio A and for 
Portfolio B, respectively. In practice, it is recognised that risk measures on the aggregated risk are sensitive to 
the strength of the dependence between individual risks. Actually, by taking into account the dependence within 
and accross the ceding insurer portfolios which is determined by the parameters in Table \3.S\. the aggregated 
risk i ?2 of the reinsurer is riskier than in the independence case. Therefore, based on VaR and TVaR as a 
risk measure, the reinsurer needs much more risk capital in the dependence case. Furthermore, for a different 
confidence level p, it can be seen that the deviation from the independence assumption is greater for VaR than 
for TVaR. 


Confidence level 

Independence case 

Laplace case 

FCM case 

V 

VaRn^ip) 

TVqRr^ (p) 

VaRR^{p) 

TVoRr^ (p) 

VaRR^ip) 

TVaRR^ip) 

90.00 % 

11.73 

22.64 

11.98 

22.93 

13.92 

25.35 

92.50 % 

14.98 

25.76 

15.24 

26.06 

17.36 

28.62 

95.00 % 

19.47 

30.10 

19.75 

30.41 

22.10 

33.14 

97.50 % 

26.97 

37.40 

27.27 

37.73 

29.93 

40.68 

99.00 % 

36.64 

46.85 

36.97 

47.21 

39.93 

50.40 

99.50 % 

43.80 

53.89 

44-15 

54.25 

47.30 

57.59 

99.90 % 

60.08 

69.92 

60.45 

70.31 

63.91 

73.89 


Table 4.1: Deviation of VaR and TVaR from the independence case. 


It is well known that risk diversification across portfolios arises from aggregating their individual risks, see 
e.g., HEi. In this respect, by considering the TVaR as a measure for the risk capital, the diversification 

benefits Dp(Ti,T2) are quantified as the relative reduction of the risk capital required for the whole portfolio of 
the reinsurer from aggregating the stop loss risk Ti and T2 as follows 

(rj, _ TVaRR^jp) 

TVaRTAp) + TVaRT,{p)' 

As presented in Table 14 diversification benefits increase with the confidence level. Conversely, the deviation 
from the independence case yields a reduction of the diversification benefits which is obvious since the full 
diversification effects are attained when risks are independent. 



P 

TVaRR^ip) 

TVaRrM 

TVaRT2 (p) 

Dp(ri,T2) (%) 

Independence case 

95.00 % 

30.10 

24.87 

18.26 

30.19 

97.50 % 

37.40 

32.34 

24.26 

33.92 

99.00 % 

46.85 

41.89 

31.97 

36.56 

99.90 % 

69.92 

64.84 

50.55 

39.40 

Laplace case 

95.00 % 

30.41 

25.11 

I8.4I 

30.13 

97.50 % 

37.73 

32.58 

24.42 

33.82 

99.00 % 

47.21 

42.14 

32.13 

36.44 

99.90 % 

70.31 

65.07 

50.72 

39.28 

FCM case 

95.00 % 

33.14 

27.71 

18.64 

28.52 

97.50 % 

40.68 

35.40 

24.69 

32.29 

99.00 % 

50.40 

45.14 

32.44 

35.03 

99.90 % 

73.89 

68.32 

51.08 

38.11 


Table 4.2: Diversification benefits based on TVaR of the aggregate risk R2 and the individual risk Ti,i = 1, 2. 











































4.2 TVaR capital allocation 

In this section, we derive analytical expressions for the amount of capital allocated to each individual risk of the 
reinsurer under the TVaR principle. In the enterprise risk management framework, to absorb large unexpected 
losses, reinsurer are required to hold a certain amount of economic capital for the entire portfolio. In this 
respect, the so-called capital allocation consists in attributing the required capital to each individual line. This 
allows the reinsurance company to identify and to monitor efficiently their risks. In the literature, many capital 
allocation techniques have been developed, see for instance IT^ [Ml l4l IT^ and references therein. In practice, 
it is well known that the TVaR principle takes into account the dependence structure between risks and satisfy 
the full allocation principle. More precisely, if Ti is the aggregate risk where Ti is a rv with finite 

mean that represents the individual risk of the reinsurer, the amount of capital TVaRp(Ti,Rn) required for 
each risk for i = 1,... ,n, is defined as 


TVaRp{Ti, Rn) 


1-p 


where p G (0,1) is the tolerence level. The full allocation principle implies 


TVaRjiJp) = J^TVaRp(T,,R„) 


(4.3) 


which means that, based on TVaR as a risk measure, the capital required for the entire portfolio is equal to 
the sum of the required capital of each risk within the portfolio. The following proposition develops an explicit 
form for TVaRp{Ti, R 2 ), * = 1,2, in the case of stop loss mixed Erlang type risks. In addition, we define below 
Si.k, S2,k, S2,k as in Proposition 14.II and we denote Xp := VaRa^ip). 


Proposition 4.4 Let (Vi,..., X2k) SME{/3, Q) with (32k > A, * = 1,..., 2fc — 1 and dj > 0, j = 1, 2. If 
further Tj,j = 1,2 has finite mean then 


1 


TVaRp{Ti,R2) — — 


P + Us^^k(xp,di,d2,Z{P2k)) 

2k I . 

l3„^ ( Ass,;. (^2)17 {Xp, di,Z{P2k)) 


1=1 ji 


m=l 


+u g (xp, di, d2, Z{P2 k))) 


Example 4.5 In this example, we consider the same individual risks and dependence parameters as in Example 
\3.}^ and the reinsurance programs as in Example \4.3\ Based on TVaR as a risk measure for quantifying the 
risk capital required for the whole portfolio the required capital of each stop loss risk Ti,i = 1,2 are evaluated 
for different confidence level p. Since Ti is riskier than T2, as displayed in Table more eapital is required 
for Ti eompared to the amount needed for T2 ■ 
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Laplace case 

FGM case 

V 

TVaRn^ip) 

TV aRp [Ti , i?2) 

TV aRp {T 2 , i?2) 

TVaRn^ip) 

TV aRp (Ti , i?2) 

TV aRp {T 2 , i?2) 

90.00 % 

22.93 

14.56 

8.37 

25.35 

16.19 

9.16 

92.50 % 

26.06 

16.85 

9.21 

28.62 

18.62 

10.00 

95.00 % 

30.41 

20.15 

10.26 

33.14 

22.12 

11.02 

97.50 % 

37.73 

25.99 

11.74 

40.68 

28.21 

12.47 

99.00 % 

47.21 

33.94 

13.27 

50.40 

36.39 

14.01 

99.50 % 

54.25 

40.03 

14.22 

57.59 

42.59 

15.00 

99.90 % 

70.31 

54.20 

16.11 

73.89 

56.79 

17.10 


Table 4.3: TVaR and allocated capital to each stop loss risk = 1,2, under the TVaR capital allocation 
principle. 


4.3 Reinsurer default analysis 

In the enterprise risk management framework, reinsurers are obliged to hold a certain amount of capital K > 0 
in order to be covered from unexpected large losses. The amount of this capital is determined so that the 
reinsurer will be able to honor its liabilities even in the worst case with high probability. For instance, in the 
SST, K is quantified as the TVaR at a tolerance level 99% of the aggregated risk i?„ = where Ti 

represents the individual risk of the reinsurer. This means that for 99% probability the reinsurer has enough 
buffer to pay its obligations. However, in case > K the reinsurer is in default and thus ceding insurers 
are not protected from losses exceeding K i.e. — K. By analogy to the case between the insurer and the 
policyholders, see |15j . the quantity (i?„ — K)^ is called the default option of the reinsurer or in other words the 
ceding insurers deficit with U{K) := E{(R„ — Rf)+} the value of the default option. In view of the full capital 
allocation principle, for a given risk capital K required for the entire portfolio of the reinsurer, if Ki, i = 1,..., n 
is the risk capital needed for each individual risk then K = Furthermore, the value of the default 

option is also defined as the sum of the value of the unpaid losses U{Ki,K) := E {(Ti — of each 

ceding insurer(s) reinsured lines of business, specifically (see e.g. [4]) 

n 

U{K) = Y,U{K,,K) 

i=l 

Next, we determine U(K) and also derive an explicit formula for the default probability of the reinsurer 4>{K) := 
P(i ?2 > K) where i ?2 = Ti + T2. Furthermore, given that the reinsurer is in default, analytical expressions of 
the unpaid excess losses of each line of business of the ceding insurer(s) are derived. 


Proposition 4.6 If [Xi ,..., X 2 k) SME{f3, Q) with P 2 k > A, * = 1, • ■ ■, 2fc — 1 and dj > 0,j = 1,2 then for 
a given risk eapital K > 0 

U{K) = {diWs,,, {K, d2, Z{p2k)) + {d2)Us,,, {K, di, Z{P2k)) 

\ 2k I / 


FUsi^k+S2,kiFiTdi,d2, Z{fS2k))j + ^ ^ ^ ^ 7im 

'' i=l . j'i m=l 

+Ps2.k (^ 2 )%,, {K, di, Z{f)2k)) + (K, dud2, Z{f32k))^ 


[Fsk^d^)Us2^K,d2.Z{hk)) 

-KFn,{K), 


where F^^{K) = 1 — F]i^{K) with Ffl 2 (.) is defined in Proposition\4-l\ 
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Remarks 4-7 In view of Proposition ED analytical expression for the default probability of the reinsurer is 
given by 


= 1-^1 yFs, i^{di)Fs^ ,^id2 +K) + „(di + s)Fs^^i^{d2) + (di,^2,FT) 

2k I , 

+ W_"ioTn\Fg^j^{di)Fg^^{d2+K) 

1=1 3l,h,---,3l ™=i ^ 


(4.4) 


Proposition 4.8 Let Ki,i = 1,2 6e i/ie capital required for each stop loss reinsurance portfolio of the reinsurer 
such that K = Ki + K2. Given that the reinsurer is in default, if (Xi,... ,X 2 k) ~ SME{(3,Q) with fi2k > 
Pi,i = 1,... ,2k — 1, dj > 0, j = 1, 2 and Tj,j = 1, 2 has finite mean, then 

U{K,,K) = ^,(^Fs,Jd2)Us,jK,di,Z{P2k)) + Us,AK,dud2,Z{f32k))^ 

2k I . 

1=1 31 , 32,■■■,31 in=l ^ 

+Ug^ jK,dud2,Z(/32k))^ -KiFrAK). 


Example 4.9 Consider the required capital for the entire portfolio K as TVaRR2(p) o,nd for the individual 
risk Ki as TVaRp{Ti, R2),i = 1,2 presented in Table From Table E3 we can see that an increase of the 
confidence level p yields an uprise of the capital required for the entire portfolio of the reinsurer which in turn 
decreases the default probability and also the value of the unpaid losses of each portfolio of the ceding insurer. 



V 

K 

m) 

U{K) 

Ki 

U{Ki,K) 

K2 

U{K2, K) 

Independence case 

95.00 % 

30.10 

0.01860 

0.19288 

19.69 

0.15436 

10.41 

0.03852 

97.50 % 

37.40 

0.00929 

0.09483 

25.47 

0.07928 

11.93 

0.01555 

99.00 % 

46.85 

0.00370 

0.03725 

33.35 

0.03228 

13.50 

0.00497 

99.90 % 

69.92 

0.00036 

0.00360 

53.59 

0.00321 

16.33 

0.00039 

Laplace case 

95.00 % 

30.41 

0.01863 

0.19338 

20.15 

0.15586 

10.26 

0.03752 

97.50 % 

37.73 

0.00930 

0.09750 

25.99 

0.07979 

11.73 

0.01525 

99.00 % 

47.21 

0.00371 

0.03731 

33.94 

0.03237 

13.27 

0.00494 

99.90 % 

70.31 

0.00037 

0.00360 

54.20 

0.00320 

16.11 

0.00040 

FGM case 

95.00 % 

33.14 

0.01870 

0.19924 

22.11 

0.16237 

11.02 

0.03687 

97.50 % 

40.68 

0.00932 

0.09740 

28.21 

0.08212 

12.47 

0.01528 

99.00 % 

50.40 

0.00372 

0.03805 

36.39 

0.03286 

14.01 

0.00519 

99.90 % 

73.89 

0.00037 

0.00364 

56.80 

0.00317 

17.10 

0.00047 


Table 4.4: Default probability, default value option of the reinsurer and unpaid losses of the insurer. 




























11 


5 Proofs 


Proof of Proposition [37l] The joint tail probability of {Si.k, S 2 ,k) is determined in terms of the joint density 
of {Xi,..., X 2 k) as follows 


P(S'i,fc > Ml,5'2,fe > U 2 ) 




h{x)dxi ... dx 2 k- 


(5.1) 


Refering to dm, the joint density of (Xi,..., X 2 k) is given by (set 7 ^-,^ := E {gj^ (^j;.)}) 

2k y 2k h 

Hx) = n 1+XI X n ) 

i=l h=2 l^ji<j2<...<jh^2k i=l 

2k 

= fji{xjj^) /i(^i) "b , 12 /jl (^Jl )/i 2 (^l 2 ) fi{Xi) 

»=i jeCXlii} ieC'\{jij2} 

2k-l 2k 

+ . . . + fjiiXji)fj2kiXj2k) + ^ji,...,j2k 

i=l i=l 

where fiixt) = g{xi)fi{xi), 

2k 

Cl = 1 + X X 7ii7l2 - X X X “ii.i2.i3 7n7j2 7t3 + ■ • ■ + (-l)^''ai,...,2fc H 


31 32 33 


Cjl ~ X! ( X! “ 11 .12 712 + X! X! Q^ll.l2,j3 7l2 7l3 +■■■ + ( 1) 0-1,...,2k 


11 12 


ieCMii} 


C11.12 = XX( “11.12 “X '^11.12 .13 7i 3 + XX “11.12.13 .I4 7i3 7i4 + ■ • ■ + (-l)^^'ai.....2fc 

11 12 ^ 13 13 14 ieC\{jlj 2 } 

1.12.13 “ X! “11 .12 .13 .14 7l4 + X! X! “11 .12 .13 .14,15 7l4 7l5 + ■ • ■ + (~1) Q!l,...,2fc 7i 

14 14 15 iGC\{jl,l2.l3} 


Cii.i2.i3 “XXX 


,71 


Cll.....j2fc-1 ^ ^ ^ ^ ‘^11 .■■■.l2fc-l '^1.....2fe7l2fc I 

11 l2fc-l 

Cjl.....l2fc “ Q^l,...,2fci 


with C = . .,2fe}, ji e C',J 2 € C'\{ji},j 3 G C'\{ji, ^ 2 }, • ■ •, J 2 fc G C'\{ji,... ,j 2 fc-i}- 

After some rearrangements, one can express h{x) as follows 


2k 2k I 

Hx ) = Cin.^*(^*)+x X Cii.i2.....ii n X(2;ij n 

i=l i=l 11.12.■•■.ji h=l it{h,j2,...,jl} 

2k 2k 2k 

~ Cl TT fijxi) + X! X! Cji. 12 .•••.ji /i ( 2 ^ 1)1 

i=l 1=2 hd2,...,3l i=l 


where for i = 1,. ■., 2fc 




Mxi) if * ^ {,7’i,j2,---,ji}, 

f,{xi) if i G {ii,j2, ■ • ■, ji}- 


By Lemma IA.41 fi is a pdf of a mixed Erlang distribution, therefore one can write (I5.1I) as a sum product of 
convolutions of mixed Erlang risks as follows 


F{Si,k > ui,S'2,fc > U 2 ) 




2 


k-1 

fi{xi)Fk{ui -xi- Xk-i)dxk-i ■ 


.. dxi 
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QO pOO 


2k-l 


■■■ H Mxi)Pk{u2 - Xk+l - ■ ■ ■ - X2k-l)dX2k-l ■ ■ -dXk+l 

/U2 ^U2—Xk + 1 '^U2 — Xk + l—--- — X2k-2j^—}^j^l 

2k I 

7im 

i=i iij2.---.il >Ti=i 

p(X) pco poo k—1 

X / ... Y\_ fiixi)Pk{'‘J-i - Xi - ■ ■ ■ - Xk-i)dxk-i ■.-dxi 

Jui Jui-Xi -'Ui-Xi — ...-Xk- 2 i^l 


poo poo 


2k-l 


fi{xi)Fl{u 2 - Xk+i - . ■. - X 2 k-i)dx 2 k-i ■ ■-dxk+i- (5.2) 

'U2 ^U2—Xk + 1 ^ U2 — Xk + l—--- — X2k-2 

Provided that /32k > /3i,i = 1, ■■■ ,2k — 1, hy Lemma lA.51 each *—th mixed Erlang component of (15.21) can be 
transformed into a new mixed Erlang distribution with a common scale parameter Z{j32k)- Therefore, with 
the help of Remark IA.7[ by convolution (15.21) can be expressed as a sum product of two mixed Erlang survival 
function as follows 

2k i 

f{Si^k>Ui,S2,k>U2) = 6^Si,fc(Ml)T’S2,ic(w2) + ^ n 

ii.i2.....ii m=l 

Thus the proof is complete. □ 

Proof of Proposition 14.11 Similarly to the independence case described in Lemma IA.21 the df of i ?2 is of 
mixed distribution and can be expressed in terms of the joint df of {Ti,T 2 ) as follows 


Pr2 (s) = 


P(ri=0,r2 = 0) for s = 0 

P (Ti = 0,0 < T 2 < s) + P (0 < Ti < s, r 2 = 0) + P (Ti + T 2 < s, 0 < Ti ^ s, 0 < r 2 ^ s) for s > 0 

f Psi,k,S 2 ,kidi,d 2 ) for s = 0 

1 Fsi,k,S 2 ,kidi,d 2 + s) + Fs^^^,S 2 ,kidi + s,d 2 ) + + di,s + d 2 ) for s > 0. 

By Proposition 13. II and Lemma FA. 2 1 Fr^Fs) can be written in two terms as follows: 

• the discrete term 

2k i 

Psi,k,S2,kidi,d2) = ^iFs^^k{di)Fs^^^{d2)+ II(<^ 2 ), 


3l,32,...,3l 


m—1 


• the continuous term 

, 82 ,k (*^1 ^2 + s) 


2k I 

= 6 ^'Si.,(dl)CS 2 ,fc(c «2 + s) + X] n + s) 

jij2,---di rn^l 

2k I 

+^iPsi,kidi + s)Fs2^k{d2) + y~! y~! ^3i,32,...,3i n 'yjm 


1^1 jl,j2:---Jl 


2k 


HlFs^,,.,+S2,Adl,d2,s) + J2 ^n,32.-,3i ]l"/jmFs,^^+s2,^{di,d2,s). 


1=1 31 , 32 , ■■■,31 


This completes the proof. 

Proof of Proposition [431 In view of (14.31) 


TVaRp (Ti, i?2) 


^{Tll{R2>VaRR,,{p)}) 

1-p 


1 P JvaRn^ip) 


E(Til{R2=4)dS- 


□ 


(5.3) 
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First, we need to calculate E(ril{fl 2 =s}) as follows 

poo 

E(ril{^2-s}) = / ufTi,Ti-\-T2^s{u)du. 

Jo 


Let 


/5'l,fc + S’2,fc (^1 5 ^2 ? ■ ^^-fSl,fe+52,fe (^1: ^2; '^)i 


fs,,,+S,,Mi^d2,u) := ■^Fs,,,+S,,M^^d2,u). 

As in Proposition 14.11 one can express E(ril{/{ 2 =s}) as follows 

^{TiI{r^=s}) = ^i(^Fs2.dd2) J ufs^^{di+u)du +J ufs^ ^+s^^{di,d2,u)du^ 


2k I 

Y[l3m(Fs^^{d2) l ufs^^{di+u)di. 

1=1 ii j2,... j'l 


m—1 


+ X '^fsi,k+S2,M^’d2,u)duj. (5.4) 

By Lemma lATTl for Xi ^ ME{P, Qi) and di > 0, i = 1, 2 

pS 1 

/ u/x(dz + u)dw = - + l)Afc(dj,,5, Qi)lFfc+2(s,/3) =: [/xi(s,di,/?). (5.5) 

•^0 P k=0 

Similarly, by Lemma I A. 2 1 

rs OO OO 

/ ufxi+X2idi,d2,u)du = — + 0^k{di, l3,Qi)Aj{d2, l3,Q2)Wk+j+3{s, f3) —: Uxi{s,di,d2, P). (5.6) 

do P k=0 j=0 

Taking (15.51) and (15.611 into account, one may write (15.4p as follows 

]E(T'il|fl.2=s}) = ^i[^F’s2_^(d2)t/si,fc(s,di,2'(/32fc)) + Usi_k+S2,k{s,di,d2, Z{(]2k)) 


2 k I 

F y~! y~! It Ijrr. ( ^S2,fc (d2)U (s, di,Z{P2k)) 

1=1 31 , 32 ,■■■jl rn=l 


FUSiS^^dl,d2,Z{P2k)) ]• 


(5.7) 


Therefore, refering to (15.3p (set Xp := VaRp{R2)) 


TVaRp{Ti, R2) — 


1 


1 -p 
2k 


ii[Fs2,k{d2) Usi,ds,di,Z{P2k))ds + Usi^kis,di,d2,Z{l32k))ds 


Zfi, L / ^CXD 

+ X] E ^3i,32,---,3l n '^J™(^S2,fc(^2) / U§^^^{s,di,Z{l32k))ds 

1=1 31 , 32,■■■,31 ni=l ^ “p 

l-OO _ ^ 

+ / Ug^^^{s,di,d2,Z{P2k))ds 


Hence, the result follows easily. 

Proof of Proposition 1X61 Bv definition 


□ 


UiK) = E {(i?2 - K)+} = Fr, iK)E {R2 - K\R2 > K} . 

Hence in view of Remark 14.21 

U (K) = ^i( {di)Us2,k {K, d2, Z {132k)) + Fs^^^ {d2)U {K, di,Z (^2fc)) + Usi,k+S2,k {K, di,d2,Z (/32k))) 
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2k I , 

1=1 jl,j 2 ,---,jl rn=l ^ 


+Fs.Jd2)Us^^^{K, di,Z{f32k)) + (K, di,d2, ZiP^k)) j - KFr,{K), 

establishing the proof. 

Proof of Proposition [478] The unpaid losses of the ceding insurer line of business is defined as follows 

roc 

U{Ki,K)=E{{T,-Ki)^r^^k}} = / E{Ta{R,=s})ds-K{FR,{K). 

Jk 


In light of (15.71) 


U{K,,K) = ^AFs,Jd2)UsUK,duZ{p2k)) + UsUK,di,d2,Zip2k)) 


2k I , 

+ E E J2,...,ii 11 f ^Sa.fc ^ (if, di, Z'(/32fc)) 

1=1 h,n,---,3l m=l ^ 


1=1 31,32,■■■,31 

+% ^ (if, di, d2, Z[P2k ) ) ) - (K) . 


Hence the proof is complete. 


□ 


□ 


Appendices 


Appendix A Properties of mixed Erlang distribution 

Lemma A.l For a deductible d > 0, if X Q) then the df ofY := {X — d)_|_ is given by 


Fyiv) = 


Fx{d) for ?/ = 0, 


where 


with 


Fxiy + d) for ?/ > 0, 

oo 

Fx{y + d) = ^ Afc(d,/3,Q)ITfc+i(2/,/3), 




^ oo 

Ak{d,P,Q) = - '^qj+k+iWj+i{d, P). 

^ 3=0 


(A.1) 


(A.2) 


Lemma A.2 Let Xi and X 2 be two independent risks such that Xi ^ ME{/3,Qi),i = 1,2. If di,i = 1,2 are 
positive then R 2 = Yi + Y 2 , with Y^ = {Xi — di)+,i = 1, 2, has df 

/or , = 0. 

I Exi{di)Fx 2 {s + d 2 ) + Fx 2 {d 2 )Exi{s + di) + Exi+X 2 idi,d 2 ,s) for s > 0, 


where 


Fxi+X 2 (di, ^ 2 , s) = EEa fe(di, /3, Qi)Aj{d2, (d, Q2)Wk+j+2(s, P). 

k—0 i=o 


15 


Remarks A. 3 Given the tractable of the df of R 2 , the VaR of R 2 at a confidence level of p € (0,1) is the 
solution of 

Fxi {d{)Fx^{d2) + Fxi {di)Fx:^{VaRR^{p) + ^ 2 ) + Fx^{d2)Fx^ {VaRR^{p) + di) + Ft^ {VaRn^ip)) = P, 
which can he solved numerically. In addition, the TVaR of R 2 at a confidence level p € (0,1) is given by (set 

Xp := VaRR^{p)) 

13- 


TVaRii^{p) = {Fxjid.i)U Xi{xp,di, P) + Fx2{d2)U xAxp,d2, (3) + U Xi+x-i[xp,di,d2, (3) 


where 


Uxiixp,d,^) = ^(fe + l)Ak{di, l3,Q,)Wk+2{xp, 13), 




00 00 


Uxi+X2{xp,di,d2, (3) = + j + 2)Afc((ii, /3, Qi)Aj{d2,/3, Q2)Wk+j+3{xp, /3). 

k—0j—0 

Proof. Since Yi and Y 2 are independent risks which have mixed distribution, the df of R 2 can also be expressed 
as a df of a mixed distribution which depends on the value of s as follows: 

• the discrete part of is obtained for s = 0, specifically we have 

Fr, (0) = P (Yi + F 2 s: 0) = P (Yi + Y 2 = 0) = P (Yi = 0, Y 2 = 0) = Fx, (di)Fx. (^ 2 ), (A.4) 


• for s > 0 the continious part of Fr^ is given by 
Fr,{s) = P(Yi+r2^s) 

= P (Yi + F 2 ^ s, Fi = 0,0 < Ya < s) + P (Ai + F 2 ^ s, 0 < Fi ^ s, Y 2 = 0) 

+P (Yi + Y 2 ^ s, 0 < Yi ^ s, 0 < Y 2 ^ s) 

= P (Yi = 0,0 < Y 2 ^ s) + P (0 < Yi < s, Y 2 = 0) + P (Yi + Y 2 ^ s, 0 < Yi ^ s, 0 < Y 2 ^ s) 

= Fxi{di)Fx2{s + d2) + Fx2{d2)FxAs + di) + Fx^s - u + di)fx2{u + d2)du. (A.5) 

Jo 

Let Ft2 (s) := Jq Fxi (s — u + di)fx2{u + d2)du, by Lemma lA.il this can be written as 

00 00 

Pt 2 (s) = Afc(di,/3i)Aj(d2,/32) / Wk+i{s - u,l3)wj+i{u,j3)du. (A.6) 

k=oj=o do 

It can be seen that lYfc+i(s — u, (3)wj+i{u, I3)du is a convolution of two independent Erlang risks with a 
common scale parameter /3, which is again an Erlang risk with shape parameter k + j + 2 and scale parameter 
(3. Thus combining (IA.4|) . (IA.5I) and (IA.6I) the claim follows easily. 

Lemma A.4 Let X ^ ME{f3,Q) with pdf f{x, (3, Q), if g is some positive function such that E{g(A)} < 00 , 

g{x)f{x,0,Q) 

then c{x, (3, Q) = — E{g(x)}~ again a pdf of mixed Erlang distribution with scale parameter Z((3) and mixing 

weights &{Q) = {0i,O2, ■ ■ ■), with 

00 

c(x, I3,Q) = '^ 0kWk{x, Z{P)), 


k=l 


where 


k-V 


Z{/3) = 2/3 and Ok = Ej=i I ■ _Ad3 for g{x) = 2F{x), 
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• Z{P) = j3 and 


0 


for k 


Ok = < , 

- ^rb+t) for k>t, 


Hfci 9j - 


rW 


for g{x) = X* with 1 € K, 


Z{I3) = (3 + t and Ok = -si with fi = for g{x) = e with t gN. 


Proof. We have 


c{x,(3,Q) = 


9{x)fix,P,Q) 


E{ffW} 

For g{x) = x* one can write (IA.7I) as follows 

c{x,j3,Q) = 






1 ^ 


^t+fe-lg-/3x 


E 


E oo 

1=1 'll 




TUT 


/ r(s) 

f 

°° „ r(j+t) 

ro) 


E 

/ ^ y Y^oo 
s=i+l ^ Z^?=l 


Wsix.P) 


'^OsWsix^P), 


with 


(A.7) 


0 


for s €t. 


0 . = < , 

"V^o for s<t. 

Zj=i 9i r(j) 


For g{x) = e (IA.7I) can be expressed as follows (set f3 := -^) 

c{x,P,Q) = 


For g{x) = 2F{x), see [2] for the proof. 

The results presented in the next two lemmas can be found in Section 2.2 of [24] and Section 7.2 of m, 
respectively. 

Lemma A.5 If X ME{(3i,Q), then for any positive constant (32 > f3i we have 

X ^ me{P 2 . $(g)), mq) = (V'l,^^ 2 ,.. o, 


1 


E 


qk- 


pk 


•rk \ 


2.fc-lg-(/3+i)a 


E 


9fe/3 


E OO n-I 

1=1 qjp 


\wk{x,l3 + t) 


'^0kWkix,/3 + t). 


k=i 


where 


k 

= E 

2 = 1 




k>l. 
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Lemma A.6 Let Xi,X2 be two independent random variables such that Xi ~ ME{j3,Qi)^i = 1,2, then S2 = 
+ A 2 ~ ME{P, n{Qi, Q 2 }) with 


■^i{Qi,Q2} 


0 for 1 = 1, 
Ei=i 9ij' 92,/-i for I > 1. 


Remarks A.7 According to Cassette et al. (2012) (Remark 2.1), the results in Lemma I A. 61 can be extended 
to Sn = provided that Xi ,..., A„ are independent, Xi ^ ME{fd, Qi) for i = 1,..., n. Specifically, 

Sn ~ ME{f3, , Qn}) where the individual mixing probabilities can be evaluated iteratively as follows 


; Qn+l) 


0 for I = 1,... ,n, 

■ ■ ^Qn} qn+i,i-j for l = n+l,n + 2,.... 


Appendix B Joint density of sums of Sarmanov random vectors 


One of the main features of the Sarmanov distribution is that its pdf can be used to derive some results in 
analytical way. For instance Vernic |22j have derived general formula for the density of the sum of several rv 
joined by the Sarmanov distribution. Below we derive the joint density of n random vectors where each vector 
consists of k elements and we denote the sum of elements within each random vector as Si := X]}=(i-i)fc+i X„i = 
1,..., n. Furthermore, we assume that the joint distribution of the overall random vectors {Xi ,..., Xnk) has 
the Sarmanov distribution with any kernel function satisfying (ESI). 


Theorem B.l The joint density of {Si ,..., Sn) is given by 

n n n 

C{ui,...,Un) = nAiK) + ^ 

i=l h=2 ltSiji<j2<...<jh^n i=l 


where 

fSiiUi) = (/(i_i)fc+i * ... * fik){Ui), 


fs^{wi) = *... * fTA){ui] 




with 




fm (j-m ) if m ^ {ji,j2,...,jh}, 

(t>xm{s)fm.{Xm)ds if m & {ji, j 2 ,. ■., jh} ,m = 1,... ,nk. 


Proof. The joint density of (S*!,..., Sn) is determined in term of the joint density of (Ai,..., Xnk) as follows 


({ui,...,Un) = J ... 


h{x)dxi...dxnk-i, 


(B.l) 


—Ui ,S2—U2,...,Sn—Un 


with. X — (xi, . . . , Xnk)-) — X\ -\- ...-h Xk-j S2 — - “h ^2k 7 — ^nk—k-\-l • “h Xnk ■ 

Refering to (EU, 


nk 


Hx) = + 


2=1 

nk 


a 4 


n (xja 


h=2 l^ji<j2<...<jh^n 


fe=l 


~ ^ ^ ^ ^ 'W 4'jkiXjk)fjkiXjk) tt fm{Xm) 

h= 2 l!^ji<j 2 <...<jh-^n k=l 'm({ji,j 2 ,...,jh} 


i=l 
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nk n nk 

?i=2 l^ji<j2< - <ifc^n m=l 


i=l 


where for m = 1,..., nfc 




fmiXm) if {ji,j2,---,jh}, 

(l>miXm)f7n{Xm) if m G {jl, j2, ■ ■ ■ , jh}- 


Therefore, one can express (IB.II) as a sum of convolutions as follows (set Xik := (x(i_i)fc+i,...,Xifc_i),* = 

1,... ,n) 




2=1 




ik — 1 ik — 1 

n Xj)dxik 

T,— {i—l)k-\-l m—{i — l)k-\-l 


/i=2 l^ji<j2<---<jh^n 

Up p ik—1 trx. — ± 

n/ / n f^nKx„^)fl^\uj - Xm)dx 

1—1 m=( 2 —1)A;+1 m =(2 —1)«;+1 

(*), 


ik—1 


= + aii.-.j;. K)- 


i=l 


h=2 i^ii<i2<...<ih^" 


i=l 
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